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I.   OBJECTIVES  AND  SUMMARY 

This  dissertation  treats  the  examination  of  the  effect 
on  fluid  flow  of  the  theory  of  non-Newtonian  fluids.  The  gen- 
eral sethods  of  Reiner  and  Rivlin,  as  developed  by  Truesdell, 
are  followed. 

An  investigation  is  made  in  Chapter  4.  of  the  func- 
tional dependence  of  the  pressure  on  the  space  variables  and 
its  effect  on  the  possibility  of  existence  of  a  particular  flow. 
This  functional  dependence  is  closely  related  to  the  form  of  the 
second  order  viscosity  coefficients.   In  many  cases,  making  the 

sane  assumption  of  such  dependency  of  pressure  as  is  made  in 

> 
the  linear  case,  together  with  taking  the  viscosity  coefficients 


as  constants,  forces  only  the  linear  solution.   This  may 
lead  to  the  erroneous  conclusion  that  the  particular  flow  of 
the  non-Newtonian  fluid  is  not  possible.   This  functional 
.dependence  Is  discussed  in  the  general  case  In  Chapter  4.  and 
in  many  of  the  particular  flows  studied. 

The  stress  tensor  is  developed  in  such  a  manner  as  to 
emphasize  the  Kelvin  and  Poynting  effects  of  the  non-Newtonian 
theory . 

Several  specific  flows  are  considered.  These  flows 
were  chosen  as  a  result  of  many  considerations,  the  primary  one 
of  which  was,  of  course,  that  new  results  could  be  obtained. 
Other  considerations  prompting  the  selection  of  specific  flows 
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were:  (1)  the  form  of  the  rate  of  deformation  tensor,  (2)  the 
coordinate  system  in  which  the  flow  could  be  considered,  (3) 
variation  of  observable  phenomenon,  and  (4)  solvability  of 
the  equations  of  motion. 

Chapter  2.  reviews  the  symbols,  definitions  and  equations 
which  are  further  developed  in  the  Appendix.   The  importance  of 
complete  understanding  of  physical  components  of  tensors  and 
their  relationships  to  the  other  tensor  components  cannot  be 
overemphasized  when  computations  are  necessary. 

Chapter  3.  reviews  the  definitions  of  fluids  and  develops 
the  general  form  of  the  extra  stress  tensor  which  is  used  pri- 
marily herein  for  the  explicit  display  of  the  terms  contributing 
to  the  Kelvin  and  Poynting  effects. 

Chapter  4.  discusses  the  latter  effects  and  the  general 
effect  of  the  pressure  on  the  non-linear  theory. 

In  Chapter  5. ,  simple  Poiseuille  flow  is  studied.   By 
means  of  inversion  of  sories ,  a  general  solution  is  provided. 
An  explicit  form  of  the  pressure  is  also  produced.   It  is  shown 
that  the  pressure  is  dependent  on  the  downstream  variable,  a 
phenomenon  not  observed  in  the  linear  case. 

Poiseuille  flow  of  a  non-Newtonian  fluid  in  a  circular 
pipe  is  completely  solved  in  Chapter  6.   In  addition,  a 
generalization  of  the  Hagen-Poiseuille  Efflux  Law  is  given. 
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I  in  Chapter  7. ,  we  demonstrate  the  existence  and  study 
the  structure  of  shock  waves  in  non-linear  continuum  theory.  The 
steady  one-dimensional  flow  of  a  trivariate,  heat  conducting 
fluid  with  non-linear  viscosity  is  considered  in  the  absence 
of  external  forces.   The  discussion  is  limited  to  perfect  gases. 
However,  the  proof  of  existence  carries  over  to  fluids  with 
more  general  equations  of  state.   Geometrical  properties  of 
solutions  to  systems  of  differential  equations  are  employed 
to  show  existence.   Comparison  is  made  of  thlcknessess  of  shock 
waves  in  non-linear  fluids  and  linear  fluids  after  choosing 
an  arbitrary  form  of  the  viscous  stress. 

The  steady  flow  of  a  non-Newtonian  fluid  in  a  wedge  or 
diverging  channel  is  studied  in  Chapter  8.   This  flow  provides 
an  excellent  example  of  the  computational  methods  necessary  in 
the  non-linear  theory.   We  are  lead  to  the  same  elliptic  in- 
tegral solution  as  in  the  classical  case.   The  pressure  dis- 
tribution is,  of  course,  different. 

The  problem  of  Boussinesq  is  considered  in  several  general 
cases.  The  effects  of  pressure  dependence  are  very  noticeable  in 
this  problem.  Consideration  of  the  dependency  of  pressure  on  the 
distance  from  the  boundary  is  mandatory. 

The  equations  of  motion  of  the  draining  of  a  vertical 
plate  are  solved  in  Chapter  10.   Quantitative  information  is 
provided  and  a  set  of  curves  showing  comparisons  of  thicknesses 
of  the  linear  and  non-linear  theories  with  observed  results  are 
provided. 


II.   SYMBOLS,  DEFINITIONS  AND  EQUATIONS 

General .  -  Insofar  as  is  possible,  standard  symbols  of  modern 
fluid  dynamics  will  be  used  throughout  the  text.   The  symbols 
and  conventions  of  tensor  analysis  are  used  wherever  general 
terms  or  equations  are  written.   Bach  symbol  is  defined  when 
first  introduced. 

In  actual  computations,  physical  components  are  used. 
The  beauty  and  elegance  of  the  invariant  form  of  equations 
unfortunately  do  not  generally  persist  when  the  physical 
components  are  substituted  for  the  flow  variables.  On  the 
other  hand,  physical  interpretation  of  tensor  components  is  not 
always  clear.   The  equations  using  physical  components  and  the 
comparable  ones  employing  the  invariant  forms  can  indeed  be 

quite  different  in  any  but  a  rectangular  cartesian  coordinate 

2 

system.   In  fact,  it  is  this  difference  that  lead  Truesdell 

to  write  "in  1944,  I  first  had  the  puzzling  experience  of 
getting  an  apparently  wrong  form  of  the  equations  of  elasticity 
in  polar  coordinates  by  specializing  the  perfectly  correct 

general  tensor  equation  to  the  case  at  hand.'* 

-  - -     ■        -  ■ 

C.   E.  Weatherburn,  Riemannlan  Geometry  and  the 
Tensor  Calculus.  (Cambridge":  University  Press,  T958) . 

2 

C.   Truesdell,  'The  physical  components  of  vectors  and 

tensors,"  Zeitschrlft  flir  angewandte  Mathematlk  und  Mechanik, 
XXXIII,  (1953),  345.356. = 
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We  use  bCti)   to  mean  the  physical  components  of  the 
tensor  X\   • 

We  assume,  of  course,  that  a  fluid  is  a  continuous 
medium,  with  which  the  reader  is  presumed  to  be  familiar. 
Equations  are  derived  in  the  Appendix,  but  are  given  below. 

Equation  of  Continuity. 

Let  p  be  a  fluid  density 

X    be  the  contravariant  components  of  velocity 
where    is  the  symbol  for  material  differen- 
tiation, 
then  the  equation  of  continuity  takes  the  form, 

where  "  j  j  M  is  the  symbol  for  covariant  differentiation. 


Cauchy's  Laws  of  Motion. 

Let  t    be  the  stress  tensor  (symmetric) 
f    be  an  extraneous  force  field,  and, 


then, 


X    the  components  of  acceleration, 


t% -h  d1 = e*'  <2-2) 


is  Cauchy's  first  law  of  motion. 


*See~Jlppend±x  ±or__  del ini t  iojuoX^t Q 


Energy  Equation 


Let  £  be  the  internal  energy  density,  and 
be  the  heat  flux  vector, 


then 


Defining  the  rate  of  deformation  tensor  as  usual, 
and  recalling  the  symmetry  of  U      we  obtain 

ei-t^JUi-V"  (2"3) 

the  Fourier-Kirchoff-C.  Neumann  equation. 

Extra  Stress  Tensor. 

————————————————————  , 

Defining  the  extra  stress  tensor,  Vj  ,  by 

V;  S  f»S]  +  t]  (2-4) 

where  <)*  is  the  Kronecker  delta,  and  Jp   is  any  scalar. 
Then  the  equation  of  energy  can  be  placed  in  the  form, 

Cen^lir-^-h-U-'U.-j-fi 

where  K|  is  a  local  statistical  parameter  (entropy) ,  and 

V2-L    rrz-lL  e*2!L.. 


Pv 


3n 


For  Incompressible  flow  (  g  *  O  )  ,  0\f*O  .   la  order  that  we 

may  obtain  the  sane  elegant  result  In  the  compressible  case,  we 
choose  the  scalar  IS-FT"  and  obtain, 


get)  5s  jirft^JLii  -  f*i 

a  different  forr  ri  the  energy  equation.  In  the  flows  considered 
herein,  we  assent  that  the  stress  power  is  non-negative,  that  is 

or  equivalent 2 .5  , 

Once  w>l  assume  an  equation  of  state,  say 

a  given  form, 

^ss  -  X  &»i  (Fourier's  Law) 

and  the  form  of  the  extra  stress  tensor 

VJ  *  i(  &\  >&'  f>  ~*>~>  ir,  e, ■  •  •) 


denote*?  functional  dependence  and  a,b,c...  are  constants 
of  the  fluid,  then  we  have  all  the  equations  necessary  to  study 
the  flow  of  fluids  considered  herein  and  many  others. 


4 
C.  Trues dell,  "On  the  Viscosity  of  Fluids  According 

to  the  Kinetic  Theory,"  Zeltschrift  fur  Pbysik,  CXXXI ,  (1952), 

273-289.  

|*  is  used  throughout  to  denote  functional  dependence 
rather  than  a  specific  function. 
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Classification  of  Fluids  and  Motions.  -  If  g»0  in  all  Motions 
then  the  fluid  is  said  to  be  incompress ible .   Furthermore!  if 
P     is  an  absolute  constant,  the  fluid  is  said  to  be  homogeneous 
and  incompressible « 

If  p  r.  h(0    in  all  motions!  then  the  fluid  is  said  to 
be  piezotropic. 

If  p *  #(C>*|)  in  all  motions,  then  the  fluid  is  said  to 
be  trivariato. 

If  £  *  £tt)  »  tne  motion  is  said  to  be  isochoric;  if 
-f}**-plt)>  isobar lc ;  if  £  ~  £(1*jt)  or  &  s  49 (£,  t) ,  baro tropic. 
All  other  motions  are  said  to  be  baroclinic. 


III.  FLUID  DEFINITIONS 

The  general  form  of  the  extra  streps  tensor.  -  The  definition 

of  the  extra  stress  tensor  is  a  mathematical  formulation  of 

1 
Stokes  observation  "that  the  difference  between  the  pressure 

on  a  plane  in  a  given  direction  (i.e.,  a  component  of  stress) 
passing  through  a  point  P  of  a  fluid  in  motion,  and  the  pres- 
sure which  would  exist  in  all  directions  at  P  (i.e. ,  the 
hydrodynamic  pressure) ,  if  the  fluid  were  in  a  state  of  rel- 
ative equilibrium,  depends  only  on  the  relative  motion  of  the 
fluid  immediately  about  P  ,  and  that  the  relative  motion  due 
to  motion  of  rotation  may  be  eliminated  without  effecting  the 

difference  of  the  pressure  above  mentioned."  Thus,  following 

034  5 

Reiner  ,  Rivlin    ,  and  Truesdell  ,  we  write 

t\+    tP$]    *  VJ  (3-1) 

where   v|  =0,  i}  c4j=0. 


G.  G.  Stokes,  "On  the  theories  of  the  internal  fric- 
tion of  fluids  in  motion,  and  of  the  equilibrium  and  motion 
of  elastic  solids,"  Transactions  of  the  Cambridge  Philosoph- 
ical Society,  VIII,  (1845) ,  287-319. 

o 

M.  Reiner,  "A  mathematical  theory  of  dilatancy," 

American  Journal  of  Mathematics,  LXVII,  (1945),  350-362. 

3 
R.  S.  Rivlin,  "Hydrodynamics  of  non-Newtonian  fluids," 

Nature,  CLX,  (1947),  611-613. 

4 

R.  S.  Rivlin,  "The  hydrodynamics  of  non-Newtonian 

Fluids  I,"  Proceedings  of  the  Royal  Society  of  London,  CXCIII, 
(1948),  260-ZHT: " 

5 
C.  Truesdell,  "A  new  definition  of  a  fluid,  I.   The 

Stokeslan  fluid,"  Journal  de  Mathematiques  Pures  et 

Appliquees,  XXIX,  TT95QT7'ZT5?Zfi:  " 


10 

g 

We  consider  only  isotropic  non-Maxwellian  fluids,  and 


set 


(3-3) 


U j  =  A|M  +  B&,  Ai  -+■  B*&»  JL*ha3     «-2> 

where  A  »B  and  C%  will  be  defined  below.   For  ease  in  identi- 
fying Newtonian  2nd   non-Newtonian  fluids  and  the  higher  order 
viscosity  effect,  we  re-write  (3-2)  in  a  more  convenient  form, 

where  A+jtM    is  tne  bulk  viscosity  ,  /J    is  the  shear  viscosity 
and 

K,  =  AB<a,-2/<  <3-4> 

Ka.  »  A  B^G*  . 

(3-3)  is  an  elegant  form  of  the  extra  stress  tensor 
from  which  we  shall  take  departure  in  defining  all  non-Max- 
wellian fluids.   From  (3-3),  the  higher  order  viscosity  effects 
become  immediately  evident. 

Newtonian  Fluid. 

If  all  K^sO  »  we  have  from  (3-3)  the  Newton-Cauchy- 
Poisson  form  of  the  extra  stress  tensor, 


6Ibid.   223-224. 

7 
C.  Truesdell,  "The  present  status  of  the  controversy 
regarding  the  bulk  viscosity  of  fluids,  "Proceedings  of  the 
Royal  Society,  CCXXVI ,  (1954),  59-65- 
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v]  ~  /*  3l*  5^  +  <L/4  dj 


which  Is  the  characteristic  equation  of  the  classical  theory 
of  fluids.   All  results  of  the  theory  of  linear  viscosity,  of 
course,  stem  from  this  relationship. 

8 

Stokesian  Fluid.  -  Truesdell  characterized  the  Stokeslan 

fluid  as  a  fluid  without  natural  time,  i.e.  ,  a  body  which 
offers  no  response  to  its  previous  form  or  a  "body  without 
a  memory . " 

If  we  take 

«_  Ms 
o-  p 

where  JD ,    for  incompressible  fluids  is  an  unspecified  scalar, 
while  for  compressible  fluids ,  -jb  is  the  thermodynamic  pres- 
sure, and  where  £JLn   is  the  natural  viscosity  (dimension  p=s  ), 
then  (3-4)  gives 

which  characterize  the  Stokesian  fluid. 

9 
Reiner  -  Rivlin  fluid.  -  We  continue  to  follow  Truesdell 

O 

Truesdell,  "A  new  definition  of  a  fluid,  I.,"  XXIX 
9Ibid. 


12 

by  adding  a  third  material  constant  vn     (dim  1  ) ,  obtaining  a 
fluid  with  a  natural  time.   The  |^~  characteristics  of  the 
Reiner-Rivlin  fluid  are  obtained  by  taking 

whence , 

We  observe  that  in  the  Reiner-Rivlin  fluid,  we  may 
include  \{Q    in  p    with  no  loss  of  generality.   Further,  we  ob- 
serve that  the  form  of  K|  *s  *he  same  for  both  the  Stokesian 
and  Reiner-Rivlin  fluids. 

The  form  of  the  coefficients  Gtt  *  -  We  consider  only  isotropic 
fluids  and  assume  6n  analytic.   Then  <Sy  can  be  placed  in  the  form, 

6,  -f"  a„„(Bi)*(tfiy'(Brm)" 

—  i ,  n  ■  m  •  -  -  —  —  — 

of  the  rate  of  deformation  tensor. 

The  6?x3k   are  dimensionless  functions  of  J&mt£L 

r     § 

in  the  case  of  the  Stokesian  fluid,  and  additionally  of 
**■'"    in  the  case  of  the  Reiner-Rivlin  fluid.  <fo  is  simply 
a  reference  pressure;  Q    a  reference  temperature. 
We  note  that 


1? 
\  -  A  BGoioo  (3_5) 

I 

in  both  the  Reiner-Rivlin  and  Stokesian  fluids. 


IV.  GENERAL  EFFECTS  OF  NON-LINEAR  THEORY 

The  extra  stress  tensor  in  extensio.  -  In  order  to  examine 
clearly  the  general  effects  of  the  non-linear  theory,  we  write 
the  components  of  the  extra  stress  tensor  in  extensio,  thus, 
from  (3-3) , 

•\ir|    -    A  £*  +  */A  A'i  ■+  Ko  +K.A1,  +  Kji  Ah Al 
Ml    «    A  AkH-  I/A-Al  +  K. -H<4At-f-KiA*kA* 

Vl     »     ^  Ah  +  */*•  A»   -f  K.  4- K,  A!  +  K*  &'h  &*     (4-i) 

\r  1    «  I/tAl  -nc,  A1*  ■+ K^  A'k  At 

Vi     =  2L/*  Ai  +  K,  J£  +Kt&*  A* 

tf  J    a  Z/a.  A3  +K,  &j  -hk*  Ak  A$ . 

Kelvin  Effect.  -  The  terms  K.o*"KtAt,  Ko  +  K4A«,  and  K*  *"  Ki£« 
of  (4-1)  represent  the  Kelvin  effect,  a  phenomenon  that  occurs 
in  the  theory  of  non-Newtonian  fluids  that  is  not  explained  in 
the  classical  theory.   If  a  hydrostatic  pressure  equal  to  the 
negative  of  \^9   is  wanting,  the  fluid  will  tend  to  change  its 
volume.   In  the  case  of  an  incompressible  fluid,  an  arbitrary 
hydrostatic  pressure  may  be  added  to  the  stress  and  no  new  ef- 
fect is  noted  over  the  classical  case.   This  effect  has  been 

termed  the  Kelvin  effect  by  Trueedell   after  Kelvin  who  pre- 

2 
dieted  it  theoretically  . 

1Ibid. 

2 

W.  Thomson  (Kelvin)  and  P.G.  Tait,  Treatise  on 

Natural  Philosophy  (2nd  ed.  2nd  half;  Cambr J dge :  TXnfvers i t y 
Press,  1883). 

14 
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Poynting  Effect.  -  Each  of  the  component©  \ft    »  Vx      and  yj 
of  the  extra  stress  tensor  contains  a  term  of  the  form  KiHkAt 
etc.   Hence  normal  tensions  must  be  supplied  in  the  planes  of 
motion;  if  such  tensions  are  wanting,  the  material  will  tend 
to  change  its  length.   This  phenomenon  occurs  in  the  theory  of 

non-Newtonian  fluids,  without  regard  for  compressibility  and 

3 
we  call  this  the  Poynting  effect,  following  Truesdell  . 

4 
Poynting  first  noted  this  lengthening  .   He  stated  that  "when 

a  wire  is  sufficiently  loaded  to  be  straightened,  it  is  leng- 
thened by  twisting  by  an  amount  proportional  to  the  square  of 
the  twist  and  with  a  given  number  of  turns ,  Inversely  as  the 
length."  Poynting  later  reported  the  results  of  measurements 
of  diameter  contraction,  which  verified  his  predictions. 

The  effect  of  pressure  on  the  non-linear  theory.  -  Let  us 
write  the  extra  tensor  tensor  in  the  form, 

where  the  fa     are  defined  by  the  corresponding  coefficients 
in  (3-2);  then  the  general  equations  of  motion  of  the  steady 


3Truesdell,  "A  new  definition  of  a  fluid.  I.,"  XXIX. 

4 

J.  H.  Poynting,  "On  pressure  perpindlcular  to  the 

shear  planes  in  finite  pure  stress  and  the  lengthening  of 
loaded  wires  when  twisted."  Proceedings  of  the  Royal  Society, 
LLXXII,  (1909),  546-599.  

5 
J.H.  Poynting,  "On  the  changes  in  the  dimensions  of 
a  steel  wire  when  twisted  and  the  pressure  of  distortional 
waves  in  steel."  Proceedings  of  the  Royal  Society,  LXXXV7 , 
(1912),  534-561. 
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state  flow  of  an  incompressible  fluid  may  be  written  in  the 
form 

C*i+t*i  *'rifii*f*fy  *^ *»■*>]*».       (4-2) 

Let  ft     be  constant  and  let  U,'h    and  -p    be  the  velocity  and 
pressure  which  satisfy  the  corresponding  Navier-Stokes  equation, 
that  is,  let 


We  briefly  discuss  the  conditions  under  which  U>i   also  satisfies 
(4-2) .   To  this  end  we  divide  the  pressure  in  the  latter  case 
into  two  parts,  thus  p  * f*  +  P     •   Then  66^  will  satisfy  (4-2) 
whenever 

Generally,  the  right  side  of  (4-3)  may  or  may  not  be  a  gradient. 
When  it  is,  we  can  obtain  solutions  of  this  type. 

Three  possibilities  occur  from  (4-3):  1)  f^^O.t    In  which 
case  p~p     and  the  flow  is  Newtonian;  2)  there  exists  a  non- 
trivial  and  additional  pressure  gradient,  P,»   ,  which  satisfies 
(4-3)  and  must  be  added  to  -to    in  order  for  il^    to  be  a  solution 
also  in  the  non-linear  case;  and  3)  there  does  not  exist  a 
function,  p,  ^  ,  satisfying  (4-3).   In  this  case,  secondary  flows 
may  have  to  be  introduced  in  order  to  solve  the  non- linear  flow, 
if,  Indeed,  solutions  exist  at  all  in  a  particular  case. 

It  is  further  noted  that  a  special  case  of  2)  above  occurs 
when  (AVAfj  ),£«<>• 


• 


•  ■ 
>' 


I 


*  ■ 


■''•'. 


.    ,  ■     • 
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V.   SIMPLE  POISEUILLE  FLOW 


General.  ~  This  flow  was  first  introduced  by  Euler  in  1757. 
It  is  considered  herein  to  demonstrate  the  use  of  inversion 
of  series  to  obtain  solutions. 

Let  there  be  two  infinite  parallel  planes,  the  lower 
one  fixed,  the  upper  moving  parallel  to  itself  at  a  steady 
speed,  and  assume  that  the  fluid  adheres  to  the  boundary.   We 
take  the  X  -axis  in  the  direction  of  motion  of  the  upper  plane 
and  the  UL  -axis  perpindicular  to  the  planes.   Steady  state  is 
assumed. 

Then 

» 


<A<4>  = 


o 


O 

o 


o 

o 
© 


with  elementary  symmetric  functions 

and 

The  physical  components  of  the  stress  tensor  are  given  by 


*ty)  ■  -/» 


f/^« 


"  1 

o 

a] 

o 

t 

o 

a 

o 

t 

f-y 


o   /  o 
too 

o   o  o 


too 

G    I     O 
O    O    / 


o  /  o\        f     ,*   i  o  o 
j    0O\j-   _£a4JL.  \o  I   o 

o  o  a\  +      \o  ©  o  . 
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18 
The  Kelvin  and  Poynting  effects  are  observed  even  in  the 
case  of  isochoricf low.   A  pressure  must  be  supplied  perpin- 
dicular  to  the  planes;  if  this  pressure  be  wanting,  the  upper 
plane  will  rise  or  fall,  depending  on  the  sign  of  J^    . 

The  equations  of  motion.  -  The  equations  of  motion  are 
The  solution.  -  From  (5-1) 2,  we  have 

— — 1SSJ—   —   \J    . 


hence,  from  (5-1) ^ 
where , 


Integrating  (5-2) ,  we  obtain 

the  resistance  per  area  of  the  flow. 
Again  from  (5-1) ^, 


(5-3) 


' 


V 


; 
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and  from    (5-1) 2 

Hence 

too  =  *y«;  °A~-£d'Z  (-+$'*)  <5-4> 

wherein  $#  are  constants  of  integration. 

Rivlin  considered  the  compressible  case  and,  further 
assumed  (}t f  »  Q    ,  obtaining  the  resistance  per  area 

tO*)~A  *£$'?. (-{#*) 


or 


-tox)  =  a,  %'  +  &x%'*+  a.3  <j'*7 


2 

Truesdell  has  pointed  out  that,  depending  on  bow  many 

terms  we  take  in  this  expansion,  there  will  be  an  odd  (or 
infinite)  number  of  functions  ^V  (i.e. ,  a  certain  velocity 
profile)  which  will  yield  the  same  resistance  per  area. 
However,  we  shall  assume  that  g-  is  an  analytic  function  and 
furthermore  that  tO**)    is  a  monotonic  function  of  J  .   This 
will  permit  inversion  of  the  above  series,  obtaining  q'  as 
a  function  of  U/.   We  then  integrate  term  by  term  to  obtain 
the  function  Jtlft)   •   This  will  be  accomplished  shortly. 
The  rectilinear  shearing  flow  of  a  compressible 

Rivlin,  "The  hydrodynamics ,"  CXCII1. 

2Truesdell,  "A  new  definition  of  a  fluid,  I.,"  JCXIX 
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3 
non-Newtonian  fluid  was  first  suggested  by  Reiner  and 

4  5 

Rivlin  and  elaborated  upon  by  Truesdell  .   Further  work  by 

6 
Truesdell  treats  the  insufficiency  of  the  shearing  forces 

in  the  non-linear  case.   He  also  considered  the  rectilinear 

shearing  flow  of  a  compressible  Stokesian  fluid  and  reported 

"the  possible  velocity  profiles  <fcf  now  depend  upon  the  pressure 

4p   as  well  as  the  resistance"  and  that  "as  the  pressure  is 

reduced,  the  resistance  is  diminished  as  should  be  expected 

in  a  gas.   Both  in  the  classical  theory  and  in  the  second 

approximation,  we  have  /*>*$' ?!*«*  *  &  $2*    Du*  ir*  *he  third  approx- 

imation,  the  resistance  is  diminished  by  *^ffl  from 

the  value  predicted  by  the  classical  theory  for  the  same 

velocity  profile  <V  . " 

Truesdell *s  statements  regarding  the  possible  profiles 
of  <h    are  certainly  correct.   However,  since  we  are  dealing 
with  physical  phenomena,  one  would  expect  a  more  unique  be- 
havior of  the  function  <J  . 

Returning  to  (5-3) ,  we  rewrite, 


3 
Reiner,  "A  mathematical  theory  of  dilatancy,"  LXVII. 

4 
Rivlin,  "The  hydrodynamics  >"  CXCIII . 

5 

Truesdell,  "A  new  definition  of  a  fluid.  I. ,"  XXIX. 

6 
C.  Truesdell,  "The  mechanical  foundations  of  elas- 
ticity and  fluid  dynamics,"'  The  Journal  of  Rational  Mech- 
anics &  Analysis ,  I ,  (1952) ,  125^3007 
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where  the  left  hand  side  is  a  power  series  in  a/  •  From  the 
assumptions  regarding  analytlcity  and  monotonicity »  it  is 
possible  to  Invert  (5-3)  to  obtain  J;u  t  power  series  in  «!•. 

In  either  the  Reiner-Rlvlin  or  the  Stokeslan  fluids, 
we  have 

where 


We  assume,  the    G*ojO     CODJ9tants< 
Then 


a«ft*  o/»-  x   a*  ft4  a/ 


y  s 


l^Krfl^W^W") 


and  from  (5-3)   and   (3-5) , 


A*    _  >*  a'    +■  &s£j£*f* 

-T- r3    *    an     3 


a"  A,  A,  8ft, 


+ (5-5) 


t 


■ 


••■ 


'• 
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Reverting  7   (6-5),   we  obtain 


g(|l *  _  £<*,$- A.)  -  £gi^  (A^-A.)* 


•  •  •  • 


Whereupon,  integrating,  tertt  by  tent, 

Again  we  observe  the  effects  of  assumptions  in  the 
classical  case.   For  example,  if  f\%    is  assumed  to  be  zero, 
we  immediately  obtain  a  linear  profile. 

In  (5.6),  the  condition  of  adherence  to  the  bottom 
plane  gives  A+*0;  that  of  adherence  to  the  top  plane  will 
yield  a  value  of  At    or  At.   A  further  boundary  condition  is 
necessary,  say,  Q*  at  the  upper  plane,  in  order  to  determine 
the  profile  uniquely. 

The  pressure.  -  The  pressure  was  determined  in  (5-4),  namely, 

In  the  case  of  an  incompressible  fluid,  Ke  can  be  dis- 
regarded as  previously  discussed.   In  the  case  of  a  compress- 
ible fluid,  Ko  determines  the  Kelvin  effect.   In  both  cases,  it 


7 
Quite  formally,  if, 


«  a  0,X*(fX*+C*?+<4x*  +  e>f'r 
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observed  that  the  pressure  is  linearly  dependent  on  the 
downstream  variable. 

Finally,  in  tho  case  of  isochoric  rectilinear  shear  of 


Bin       /A*)*" 

which  appear  in  the  explicit  expression  for  ^(*^)  *n  (5-6), 
indicate,  as  Truesdell  has  often  stressed,  the  appropriateness 
of  the  Stokes ian  theory  for  low  pressure  phenomena  in  gases. 
In  this  case,  there  is  no  longer  a  general  explicit  expression 
for  the  pressure  since  in  general ,  jt0  (  ?x       and  «*  are  power 
series  in  -4-   ♦   However,  in  the  case  that  fa     and  %     are 
constant , 

and  the  expression  for  the  pressure  remains  valid. 


VI.  POISEUILLB  FLOW  IN  A  CIRCULAR  PIPE 

General.  -  We  briefly  consider  laminar  flow  through  a  cylindrical 
pipe  under  steady  state  conditions  in  order  to  derive  an  extension 
of  the  classical  efflux  law.   We  take  the  axis  of  the  pipe  as  the 
2i    -axis  in  a  cylindrical  coordinate  system  /t;  fl,  z.     •  Adherence 
to  the  boundary  is  assumed.   The  velocity  field  is  given  by 

X(i)   =  (0,  o,  $10) 

and  the  acceleration  is  zero. 

The  equations  of  motion.  -  The  equations  of  motion  are  given 

The  form  of  the  pressure.  -  From  (6-1)3,  we  deduce 
dZ- 


or 


40   *  Z-k(/L)  +  k(K) 
and  therefore 

but  from  (6-1) j,  we  note  that 

it  «  f  ^/L> 
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hence  l\,=0,  rt  *  fli  »  ***<*  the  pressure  gradient  is  there- 
fore constant  downstream. 
Again  from  (6-1) 3, 

i»'fc(-*J'V--A*.  (6-2) 

Integrating  (6-1)  j_ ,  by  parts,  we  obtain 

and  substituting  for  Q'   from  (6-2) ,  we  obtain 

p  =  J.  -  A.z.  v-  -j:  r'i  +  AV*  #  ^  •  <6-3> 

The  flux.  -  The  main  flux  through  the  tube  of  an  incompressible 
fluid  is  given  by 


a=J 


X.TTfi<l<^^ 


■i 


Integrating  by  parts,  we  obtain 


a.^^C-i.^]. 


The  first  term  in  the  bracket  vanishes  because  of  the  condition 
of  adherence,  and  therefore 

a-  tn-Mflr-^- 

If  jL  is  constant,  then  from  (6-2),  we  may  write 

a' 4i* 

3      ^ 
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a  parabola,  the  classical  result  of  Poiseuille. 

Further,,  from  (6-3) 

Defining  the  inlet  pressure, 

and  the  outlet  pressure 

* 

where  X  is  the  length  of  the  pipe, 
we  have 

and,  the  flux 

Q^  ~  ___.  {-"J-*-- )  .  (6-4) 


This  expression  for  the  flux;  was  first  given  experimentally 

1  2  3 
by  Poiseuille   '  '  *      and  is  termed  the  H&gen-Poiseuille  Law  of 


1 

J.L.H.  Poiseuille,  "Experimental  research  on  the  motion 

in  tubes  of  very  small  ctiainmeter ,  *•'  academle  des  Sciences  Comet es 
Rendus,   XI,    (1840),    961-967,    1041-T04a.~ 

2 JJll.,9  »    XII,    ( 184 1 )  ,    112-115. 

3 
J.L.  tf>  Poiseuille,  "Experimental  research  on  the  motion 

in  tubes  of  very  small  diameter,"  Memoire  des  Savants  .Strangers, 

19:'  (1848)  ,  433-544. 
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Efflux  from  a  pipe.  The  formal  derivation  is  due  to  ftivlin  . 
A  generalization  to  this  law  is  derived  below  for  the  non- 
linear case. 

We  observe  that  the  velocity  distribution  and  the 
mass  flux  are  unchanged  from  the  classical  case  in  the  second 
order  (  Jv  a  constant)  theory. 


The  solution  by  inversion.  -  We  return  to  the  differential 
equa t  ion ,  ( 6-2 ) , 

±8'fc(-i8'*J--A.rt..  (6-2) 

By  methods  similar  to  those  jn  the  previous  chapter,  we 
derive 

The  mass  flux  in  the  case  of  a  non-Newtonian  incom- 
pressible fluid  is  given  by 

r  A.a*  .  /».^f».Q.*  .  .7 

s  TTtf^-—  +     4-8^  J. 

Computing  only  the  first  two  terms,  and  taking 
^-  XilXs.    ,  as  before,  we  obtain 

an  extension  of  the  Hagen-Poiseuille  Efflux  Law,  which  is 

4 
Rivlin,  "The  Hydrodynamics "  CXCI1I 
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more  conveniently  written  as 

Polseuille  used  a  "pipe"  of  radius  ranging  from  0.007 
to  0,325  nn.   The  new  law  predicts  that  if  similar  experi- 
ments were  performed  with  Hi j  Br  sufficiently  large ,  yet 
the  flow  still  laminar ,  that  a  measure  of  the  efflux  would  be 
substantially  different  from  that  given  in  the  classical 
efflux  law. 

In  general,  $(1)  will  be  an  even  power  series  in 
'*■"  ■  ■  ,  and  Q.   will  be  an  odd  power  series.  The  classical 
Hagen-Poiseuille  Efflux  Law  is  valid  whenever 

<x^e±  «  1. 


• 


- 


VII.   STRUCTURE  OF  SHOCK  WAVES 

General.  -  The  existence  and  structure  of  shock  waves  in  the 

non-linear  continuum  theory  was  established  for  the  fluids 

1 
described  herein  jointly  by  Gilbarg  and  the  author  .   In 

addition  therein,  we  computed  shock  thicknesses  in  monotonic 
gases  and  compared  these  results  with  those  of  Thomas  (con- 
tinuum theory)  and  Mott-Smith  and  Zoller4  (kinetic  theory); 
and  further  demonstrated  that  the  basic  methods  of  extending 
the  theory  to  polyatomic  gases  remains  unchanged  in  the  con- 
tinuum theory  provided  that  fluid  moduli  can  be  determined. 

That  portion  of  this  joint  work  pertaining  to  the  struc- 
ture of  shock  waves  in  the  non-linear  theory  is  included  as  it 
relates  to  the  topic  of  this  dissertation. 

The  steady  one-dimensional  flow  of  a  trivariate,  heat 
conducting  fluid  with  non-linear  viscosity  is  considered  in  the 
absence  of  external  forces.   The  flow  approaches  finite  limit 
values  at  %s.±oO    •   Such  flows  display  the  character  of  a  shock 

D.  Gilbarg  &  D.  Paolucci,  "The  structure  of  shock  waves 
in  the  continuum  theory  of  fluids,"  Journal  of  Rational  Mechanics 
and  Analysis,  II,  (1953),  617-642). 

2 
L.H.  Thomas,  "Note  on  Becker's  theory  of  the  shockfront," 

Journal  of  Chemical  Physics,  XII.  (1944),  449-452. 

3 

%H.M.  Mott-Smith,  "The  solution  of  the  Boltzmann  equation 

for  a  shock  wave,"  Physical  Review,  LXXXII,  (1951),  256-274. 

4 
K.  Zoller,  "On  the  structure  of  shock  waves",  Zeltshrlft 

fur  Physik,  CXXX,  (1951),  1-38. 
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wave  in  that  they  differ  from  their  end  states  at  X.  «■  ±  «0 

only  in  a  small  interval  of  rapid  transition, 

i 
The  velocity  field  is  given  by  X«i*  and  the  rate  of 

deformation  as  4=  -*.  U.  ,  where  the  prime  will  denote  differ- 
entiation with  respect  to  X  throughout  this  chapter. 

Equations  of  motion.  -  The  equation  of  continuity,  Cauchy's 
first  law,  the  energy  equation,  Fourier's  Law  and  the  extra 
stress,  respectively,  take  the  form 


gu.ui  *  t' 

ekt'  =  tiL'-$'  (?-1} 

f  =-**' 

XT  »  t  *f> 
Elimination  of  t    and  J-  from  (7-1)  and  integration  of  the 


resultant  system  yields 


gu? +f>  -V  =/> 


f 


(7-2) 


e*l$+t  +  fl-v"-xe'~ 


E 


where  yft,  ,  P  and  £  are  constants  of  integration  (7-2)  ex- 
presses, respectively,  the  conservation  of  mass,  momentum  and 
energy.   Eliminating  P    and  redefining 


'• 
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a. 


we  obtaira 

j^e^  is  ]z  -  ±iu,-*uf  -  c\  <7-3> 

The  form  of  the  viscous  stress*  -  We  propose  to  study  the  ef- 
fect of  non-linear  viscosity  on  the  theory  of  the  shock  profile, 
an  application  naturally  suggested  by  the  large  velocity  grad- 
ient present  in  tne  shock  front. 

Of  the  possible  generalizations  to  fluids  obeying  a  non- 
linear stress  law,  we  select  one  in  which  the  viscous  stress 
may  be  a  non-linear  function  of  the  velocity  gradient,  but 
remains  independent  of  temperature  gradients.   To  fix  our  ideas , 
we  choose  the  specific  volume, xr  ,  and  the  temperature,  0      , 
as  variables  of  state  and  define  the  viscous  stress  IT  -  V  (juJ  t  V,  O) 
with  the  following  properties: 

a.  -\f  (of  ™vri  e )  ^O  for  all  a/;  e 

b,  |L(O,vi0)  */*(v,&)>0    for  all  XT,  0 

Property  a.  asserts  the  essential  characteristic  of  fluids  that 

the  viscous  stress  vanishes  when  there  is  no  velocity  gradient. 

Property  b.  states  that  for  small  velocity  gradients,  the  stress 

is  approximately  given  by  the  linear  law  and  the  fluid  essentially 

becomes  the  classical  one. 
t  '  •  ..... 

For  convenience,  we  consider  the  function  if  ~  XT (^'j M-j O) 

which  we  identify  with  V(u!,\/'sQ)     by  means  of  (7-2) ^  and  re- 
write (7-3) 
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*  (7-4) 

For  convenience  in  exposition,  but  without  essential 
loss  of  generality,  we  limit  our  discussion  to  perfect  gases. 

However,  the  same  proof  carries  over  as  well  to  fluids  obeying 

5 
the  general  equations  of  state  considered  by  Gllbarg   in 

proofing-  existence  in  the  linear  case.   For  these  gases, 

J>-Re^,       E^CrQ,      #~cr-Cy 

where  Cm  »  £-y  »  are  the  specific  heats  per  unit  mass  at  constant 
pressure  and  volume,  respectively.   We  set  K*  ^   »  as  usual, 


then  define  %   »4(>'wj. 


Following  Becker  ,  we  introduce  dimensionless  variables 


by  means  of  the  substitutions 

/>     -   *rR©  a 


0)  =  -p-u,,   /*  »-£-->   fc>  -  -p 


and  also  define 


and  the  function 


5D.  Gilbarg,  "The  existence  and  limit  behavior  of  the 
one-dimensional  shock  layer/'  American  Journal  of  Mathematics, 
LXXII,  (1951),  256-274. 

R.  Becker  "Stossvelle  und  Detonation,"  Zeitschrift  flir 


Physik,  VIII,  (1921-22),  331-362. 
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After  the  above  substitutions,  (7-4)  becomes 

X0'  =  e  -S[(l-U))S«C]5LK§)  (7-5) 

<D  carries  over  the  properties  of  V (<**'j  i* ,  9)  •  namely 

»'  •  <?  6*/  a),  ©)  *  O  for  all  */,  £ 

b*.  ^r-L0,(JJ,h)z.  U(u)fe)  >0  for  all  *>>&. 

da)'  ' 

By  the  implicit  function  theorem,  these  properties  have  as 

an  immediate  consequence  the  existence  of  a  function  V/(ca)j§) 

such  that 

in  •some  neighborhood  of  the  parabola  AI-0,  and  such  that 
"Y»  s.0  on  this  curve.   Naturally  the  extent  of  the  domain  of 
definition  of  y>(*0*$)    cannot  be  deduced  from  the  purely  local 
properties  such  as  a*  and  b*. 

In  the  u)-§     plane  (Z. -plane)  ,  let  the  two  parabolas, 
M*0  and  L*0  intersect  in  the  points  Z,0*<fA>JW./7~f(k'o  R) 
with  u)0>U)(  .   We  see  that  such  points  could  correspond  to 
initial  and  final  states  of  a  shock  wave,  since  Z©  and  Z.t 
would  simultaneously  satisfy  (7-2)  and  therefore  satisfy  the 
Rankine-Hugoniot  shock  conditions,  namely, 
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Co  "*•  *  £>  **"  *  1H' 

tf.i-e.uz  =  //  +{.  *•*  -  /* 


e*+i«i*j?**>*-iu'''1-f 


where  the  values  of  UP   »  U-i   »  etc  are  obtained  in  the  obvious 
manner. 

Equivalent  conditions  for  existence.  -  We  will  now  reduce 

7 
(7-5)  to  a  form  wherein  the  techniques  developed  by  Gilbarg 

may  be  applied. 

We  make  the  further  assumption  that  yf  (u),§)     is  defined 

throughout  the  region ,  R,  (see  figure  1.)  bounded  by  the 

parabolas. 

Since  }£ «  O   implies  w'sO,  which  in  turn  implies  M-O 

(property  a*)  »  we  observe  that  ")£  =  0  only  on  the  parabola 

M  Lb),§) —O    .   Therefore  ^(u),  &)       is  of  one  sign  throughout  R. 

Furthermore  this  sign  is  negative  since 


Tfc 


_  ft-^s  J       -  Ik. 


&  d>   '  \  i/sO 


M'~0 


>  0 


If  we  can  prove  therefore  the  existence  of  an  integral  curve 
of  the  system 

L  (*,  S)  (7_6) 

'Gilbarg /'The  existence "  LXXII. 
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joining  the  singular  point  Z©  and  21,  ,  then  clearly  we  have 
succeeded  in  proving  it  for  the  system  (7-5). 

Proof  of  existence.  -  The  proof,  patterned  after  the  one  for 

8  9 
the  linear  case   '   rests  on  the  following  facts: 

(i)  Z0    and  7,,  are,  respectively,  node  and  saddle  point 
of  (7-6); 

(ii)  there  is  a  unique  integral  curve  SCO  of  (7-6) 
that  approaches  the  saddle,  Z-i  ,  from  within  R  as  x-*oo  . 

(lii)  S  n&s  a  negative  slope  In  R  and  cannot  intersect 
the  boundary  of  R. 

(iv)  5  cannot  terminate  in  the  interior  of  R  and  there- 
fore must  approach  Zo  • 

Steps  (iii)  and  (iv)  follow,  as  in  the  linear  case,  from 
our  knowledge  of  the  signs  of  M(t*),$)     and  L(u)j  3)      in  R. 
The  parabolas  divide  the  neighborhood  of  X\    into  four  regions 
in  which  the  slopes  of  the  direction  field  are  alternately, 
positive  and  negative,  and  (see  figure  1.)  in  R  the  slope  is 
negative,  and  in  the  direction  of  the  second  quadrant  since 
H(<J,9)  <  O    ,  LC**s6)/0     throughout  R.   We  observe  now  by  con- 
sidering the  signs  of  LCt*),§)   and  M(&j$)  that  on  the  boundary  of 
R,  all  integral  curves  are  directed  out  of  the  region  as  X  in- 
creases and  therefore  into  the  region  as  X  decreases.   Hence  if  we 
follow  the  integral  curve  StfO  backwards  from  Zi  in  the  direction 

8Ibid. 

g 
Gilbarg  and  Paolucci,  "The  structure  ,"  II 
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of  decreasing  X  (and  increasing  u)   ) ,  we  see  that  it  can  never 
leave  R  and  must  approach  the  other  singular  point  as  x— *— cO- 

It  therefore  remains  only  to  discuss  the  nature  of  the 
singular  points  at  Z.0    and  Z#  and  thereby  to  establish  (i)  and 
(ii).   The  characteristic  equation  of  (7-6)  at  Z«  and  Z-r  is 


\~%)u)      (sJa  " ~*  )**»•»•*• 


=  0 


We  observe  that 


<--«- 


a> 


*i 


similarily 


t 


2v 


while 


& 


U  fa) 


(Ml    s   ^ 


Mia(Zi) 


I*-  1 ,2. 


and 


a 


'J=\   I       -     & <z$ 


Consequently  the  characteristic  equation  is  precisely  that 
given  in  the  linear  case,  that  is, 

By  inserting  the  values  of  M    and  jL  and  subsequently 
evaluating  at  Xo    an<*  £j  »  and  after  simple  but  tedious 
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algebraic  computation,  it  is  not  difficult  to  see  that  the 
roots  of  the  characteristic  equation  are  real  and  of  opposite 
sign  at  Zj  and  real  and  of  the  same  sign  at  Zo  .  Thus  Z*  is  a 
saddle  point  and  Zo    an<*  unstable  node.   There  are  exactly  two 
integral  curves  that  approach  the  saddle  Zias  %-*<&   and  exactly 
two  that  approach  as  x~+ -  <&   •   These  pairs  correspond  to  the 
negative  and  positive  roots,  respectively  of  the  characteristic 
equation.   Two  members  of  each  pair  have  the  same  slope,  bat 
approach  from  opposite  direction.   Thus  (i)  and  (ii)  have  been 
demonstrated  and  the  existence  established. 

Numerical  computation.  -  The  numerical  calculation  of  shock 
thicknesses  in  the  linear  fluid  under  discussion  requires 
somewhat  exact  knowledge  of  the  viscous  stress ,  V (*/,<*)  0)   •  In 
an  effort  to  arrive  at  an  estimate  of  the  influence  of  non- 
linear viscosity  on  the  shock  profile,  we  consider  a  gas  for 
which  we  assume  a  quadratic  viscous  stress  of  the  form 

V*:    JL.JULU,'  J-  0C(-V,0)U,'*"  (7-7) 

where  U     is  the  usual  shear  viscosity  and  c((?S,G)   is  a  co- 
efficient still  to  be  determined.   For  lack  of  specific  empirical 
information  concerning  (K>  ,   we  arbitrarily  take  a  value  suggested 
by  the  Burnnett  approximation  of  the  kinetic  theory.   In  this 
approximation,  the  expression  for  viscous  stress  in  one-dimensional 
form  in  the  absence  of  external  forces  and  which  contribute  to  (7-7) 
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where  H>t   ,  ©t  ,  u>^    are  dimensionless  quantities  whose  values 
depend  on  the  molecular  model. 

Simply  to  serve  as  an  example,  we  take  the  values  for 
<2>,  ,  *3>a.  and  ob}    given  by  Chapman  and  Cowling 

Substituting  our  arbitrary  expressions  for  IT  into  system  (7-5) 
we  obtain 

J 


77.  77%/   —   -J - ® 


]fr  U)'   = 


—-SB— 


x9'  =   e  -5[o-w>*] 

The  positive  square  root  is  taken  for  CO     in  order  that  U)  ~  O 
imply  /f«0.   For  purpose  of  calculating  the  shock  profile,  this 
system  can  now  be  treated  as  was  that  resulting  from  the  Navier- 
Stokes  equation.   If  u,  <>*  0s    ,  then  the  expression  for  A     is 

simply  Axi*  --j"  • 

In  order  to  compare  with  previous  results  obtained  from 
the  linear  theory,  we  specialize  to  the  case  of  Helium  for  which 
S*  •44-T  and -&-*—-.   If  we  designate  by  ^/;„.   and  A^.  the 
shock  thickness  according  to  the  linear  and  non-linear  theory, 
respectively,  our  sample  calculation  show: 

afe       K#»£     &&±    -  I  *<* 

and  at   Ko"-4     ^X^l     »  »2>°  • 

^  ll-M. 

S.  Chapman  &  T.  Cowling,   The  Mathematical  Theory  of 
Non-uniform  Gftflaau   (Cambridge:  University  Press,  1939) 
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Thus  the  shock  thickness  increases  in  going  from  the  linear 
to  the  non-linear,,  the  percentage  increase  being  greater  at 
the  higher  Mach  numbers.   Of  course,  the  significance  of  these 
numerical  results  is  relative,  since  the  choice  of  the  form  of 
the  viscous  stress  was  arbitrary.   In  a  specific  real  fluid  the 
effect  of  non-linearity  might  be  greater  or  smaller  depending 
on  the  empirical  value  of  the  coefficient  of  the  quadratic 
term.   However,  the  discussion  shows  that  the  potentialities 
of  the  continuum  theory  are  by  no  means  limited  to  the  Navier- 
Stokes  Equations. 


. 


VIII.  STEADY  FLOW  BETWEEN  NON-PARALLEL  WALLS 

Introduction.  -  We  consider  the  laminar  flow  of  an  incom- 
pressible fluid  between  non-parallel  walls,  under  steady 
state  conditions.   We  take  the  intersection  of  the  walls  as 
the  Z -axis  in  a  cylindrical  polar  coordinate  system  and  as- 
sume adherence  to  the  walls.  The  coordinate, /t  ,  the  distance 
from  the  2L-axis  is  assumed  greater  than  zero.   The  equations 
of  the  walls  are  therefore  given  by  d>  *d&$, 


General.  -  Radial  flow  and  absence  of  secondary  flow  is  assumed. 


•_  ..  • 


X(H  *  (u»>o,o). 


Then  the  equation  of  continuity  takes  the  form 


5n- 


**»  — 


*c 


which  is  satisfied  by 


rt4  =%(<p,z) 


^ 


ana,  finally,  if  we  assume  that  «x—  ■  O  ,  then 


S2 


u. 


The  inertial  terms  are  therefore  given  by 
and  the  rate  of  deformation  tensor  as 

%'     *>    • 

o       o    o 
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A(iJ>  =  ■& 
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with  elementary  symmetric  functions 


The  stress  tensor  is  written 

r»  o  o*1 

%'  X*  o 
_0  o  o 

•4  ft* 

o 
o 

o 

o 
o 
o 

The  Poynting  effect  is  observed.   As  in  the  case  of  flow  in 
a  pipe,  l(ZJL)     differs  from  the  pressure  by  a  quantity  which 
is  always  greater  on  the  edges.   Hence,  if  this  pressure  is 
wanting  and  the  walls  are  terminated  at  /T.» /!,>(?>  then  the 
emerging  fluid  will  swell. 

Equations  of  motion. 


Let    ©  *  JL   . 


Then 


and  the  stress  components  may  be  written 

t(n)  -  -p  -  4?  *  £* 
#«)*  -^b  y  ^  *-  Tie 
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The  equations  of  motion  are 

_£!!«—**       i   3(**.)  ,   sew)   .    i    ay^i 

Differentiating  (8-1) x  with  respect  to  <j>  and  (8-1) 2  with 
respect  to  ft,  ,  and  subtracting,  we  obtain 

In  order  to  separate  explicitly  the  dependence  of  the 
flow  variables  on  the  space  variables,  we  introduce  several 
substitutions.   We  first  take 


Then 


£e    a  _  *e_  a*e  a   goo, 

F 

Substituting  (8-3)  and  (8-4)  in  (8-2)  and  after  some  manip- 
ulation and  gathering  of  terms,  we  obtain 


' 


• 
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-4m' »  re^'tJUHf-  ifc#+  itj(kH)'] 

+  h  i*%'  t  §."J  (8-5) 

independent  of      -fz  • 

Substituting   for    0    ,    and  rewriting,    we  have 

ifc»  H  •*<♦)  +4^HR(*>  +  £*'«*"**»'-©  (8-6) 


where 


P. W  5  [r  N*  +  I C*  H  H'  -  I  <,$  H*J 

Solution.  -  Suppose  now  that  we  take  -fti      of  the  form 

?,  =  j^a„e* 

wheue  the    <X$      are  constant.      Then 

a* 

Recalling  that      0  = A»*~    and  substituting  for    ^      and 

its  derivatives   in    (8-6) ,    we  obtain 


o=aL—  •»*• is*  *«£**•  gfr 


- 


^  *£» 


.s  + 
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or  equivalently , 


0-|;[^-)P,-*PH  +  Bri3a-^g 


Equating  the  coefficients  of  the  powers  of  1  ;rr*|  on  each 
side  of  (8-7)  we  obtain  the  following  equations, 


(8-7) 


[- 


Supposing,  in  (8-9),  that  (X?  *  O  for  all   9  >  / 


(8-8) 


(8-9) 


then 


9(^-l)P,-VP,H  +roH*-Ofor  all   ^^j 
is  a  system  of  linear  algebraic  equations  in  the  three 
unknowns p^  ,  P,  H  ,  P«  H*  • 

For  any  Q{    £    y*   ,    all  unequal ,  the  determinant  of 
coefficients  of  the  system  (8-10)  is  given  by 


(8-10) 


fi(A-t)      -p       i 
Ylr-i)      -V      i 


y*    -y    i 


the  Vanderraonde  determinant,  which  has  the  value 

and  vanishes  if  and  only  if  at  least  two  of  (£  ,  A  ,  ")/*  are 
are  equal,  contrary  to  hypothesis.  Hence  it  cannot  be  the 
case  that  all  CXy  are  non-zero. 


£*  up.  jf-*l><!i£nm 


^JW^'^i^gg1^-'^^  *"**'•  '  ■  ■•'-  iM^u-Mm*"**-**** 
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Repetition  of  the  above  argument  demonstrates  that  all  G»?  > j 
are  aero  except  for  at  most  two  different ,  say  fe(  and  fl  . 
In  this  case  (8-10)  reduces  to 


or 


where  $(  and  .#  are  to  be  determined. 

Now  from  (8-8)  ,  since  CL0  »  £  4t* 

which  is  the  differential  equation  of  the  classical  case 

12         3         4 
stuidied  by  .Jeffrey  ,  Hamel"' ,  Harrison  and  Oseen  ,  and  which 

after  integration,  yields 


•*  m  iL-fAi-  ^^  —  ^5t*  -•£*]  (8"13) 


n/' 

a      /* 

where  /h   and  Ax  are  constants  of  integration.   This  equation  is 

1 

G.8.  Jeffrey,  "The  two-dimensional  steady  motion  of 

a  viscous  fluid,"  Philosophical  Magazine,  XXIX,  (1915),  455-465. 

2 
G.  XIamel ,  "Special  motion  of  viscous  fluid,"  Deutsche 

Mathematiker  Vereinig-qng,  XXV,  (1917) ,  35-60. 

3 
W.J.  Harrison,  "The  pressure  in  a  viscous  liquid  moving 

through  a  channel  with  diverging  boundaries,"  Proceedings  of 

of  the  Cambridge  Society,  XIX,  (1920),  307-312": 

4 
„   ,      C.W.  Oseen,  "Exact  solutions  of  hydrodynaraic  dif- 
ferential equations,"  Arkiv  for  Matematik,  Stockholm,  No.  14, 
(1927);  No.  22,  (1928)": \ 


~4T 

further  studied  by  nutans  of  elliptical  Integrals. 
Ve  write,  in  lieu  of  (&-13> , 

where  two  of  the  &t,  6a.  »  &*  »  tr©  constants  of  integration 
and 

This  requires  at  least  one  of  the  &i    should  have  a  negative 
real  part.   If  we  write 

<3&l(6*}£  ®JL   (6a)>(R^(Cs) 
) 
then 

the  sign  of  equality  being  valid  only  when  the  6%  all  have 
t&e  saue  real  part,  the  different  possible  cases  have  been 
discussed  by  Hnmel. 

Returning  to  (&-11) ,  if  we  substitute  the  value  of 
from  (8-8) ,   we  have  two  equations  in  g» ,  its  derivatives  and 
their  powers.   (8-11)  is  then  a  set  of  dlaphantlne  equations 
in  (\    and  f$ .   1  have  not  been  able  to  determine  whether  non- 
sero  solutions  exist  in  this  case.  If  such  solutions  do  not 
~#*Aa^-  i1r-we*^-Hk|HE>e^ 
.possibly  including  secondary  flows ,  must  occur. 

The  pressure.  -  Suppose  that  ff  is  constant ,   i.e.^  equal  to_jfc* 
Then  by  substitution  in  equations  (8-1)  j^  £>  adding,  and   _._ 


-.-:.;•'    ;  ",'  '   :   <'  •     *  •*'      .?...■-  *'.      ' 


i  H  }     i-  ■■    ■   ■     . 


rj  .     :         '.' 


'•.    .  #*  ■■•..■•'".■■-.■• 


•    -     ■        ■  .       :.  -     •' 


• 


. 


;         ■  fc     SMI     ■•  -  •■■     '    .v  •'.-.:.  .':  ■ 

.•  •  '     •   • '       .      ••.  ''         ..  •  ■•     •• 


.• . 


»'  ' 
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collecting  terras,   we  obtain 


dtt. 


(*-*»•> +-fc|y(fr-M»  *# 


where 


We  may  therefore  assume  that 

a'fsMf— ?»©)  =^Oi(*)  (8-15) 


and 

with 

Integrating  (8-15) ^  2  with  respect  to  AC  and  <|>  ,  we 
obtain 


since  x*p-  O   from  (3-D3. 

Integrating  (8-12) ,  we  obtain  the  second  derivative 
of  the  function  CU  , 

A  f  +  +*"% + cf^/*A>  -  ° 

which,  when  combined  with  (3-14)  and  (8-1S)  yields 

Supposing  now  that  we  take,  for  convenience, 
and    Q3  3   ?/*^' 


■ 


^     v- 


* 
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then  (8-17)  become 

and  finally,  by  comparison 

_  g£fgg   ^*.  y-  fz  e  +    ••*** 

The  above  analysis  demonstrates  that  the  flow  of  a 
non-Newtonian  fluid  in  a  wedge  yields  a  velocity  profile 
identical  to  the  flow  of  a  Newtonian  fluid  provided  no 
secondary  flow  exists.   The  pressure,  however  differs  Iron 
the  classical  case  by  an  amount  f%  0   ,      It  is  noted  that  no 
restriction  has  been  placed  on  Ju  ,   while  r/  was  assumed  ccr 
stant  in  arriving  at  this  conclusion. 


IX.  THE  PROBLEM  OF  BOUSSINESQ 

General.  -  In  1330,  Boussinesq   posed  and  discussed  the 
problem  of  "how  the  retarding  influence  of  viscosity  trans- 
fers itself  throughout  the  entire  fluid  mass,  beginning  with 
the  boundary,  as  soon  as  the  stato  of  rest  ceases." 
Boussinesq  postulated  that  "at  a  given  time,  a  constant  force, 
fa,,  parallel  to  the  horizontal  axi^  is  exerted  on  a  unit  mass 
throughout  the  fluid,  in  a  manner  to  set  the  fluid  in  motion, 
except  that  the  bottom  maintains  its  own  adherence  to  the 
boundary." 

We  therefore  take  X  as  a  set  of  rectangular  cartesian 
coordinates  and  consider  an  infinite  fluid  occupying  all  space 
in  the  upper  half  plane,  X£  O  ,  and  take  the  X-axis  vertical. 
Then ,  for 

t<0;    fl=  o  ,  *;=o 

for 

X1  =  X(£)  *  (u.,  0,o) 

The  boundary  conditions  demand  that 

<j(o,t)  -  o  ,  Q(z.,o)  -  o. 

M.J.  Boussinesq,  "On  the  manner  in  which  viscosity 
enters  into  the  problem  in  a  fluid  which  goes  into  motion 
from  a  state  of  rest  and  on  its  effects  toward  preventing 
the  existence  of  a  velocity  potential,"  AcadeVie  des  Sciences 
Comptes  Rendtts,  XC,  (1380),  736-739;  9G7"imT. 
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The  inertial  terms  are  given  by 

The  rate  of  deformation  tensor  is 

r 


Sl  ivj)  * 


© 


o 


is*      © 


O 
O 


and  the  extra  stress  tensor  is 

[I       O     ©| 
O      «      0 


\o    o    o\ 
It     o    °J 

,[100]    » 3u. »,  r°  o  n     1  /ani* ,  r  1  o  «i 


with 


The  Characteristic  Kelvin  and  Poynting  effects  of  the  non- 
Newtonian  fluid  theory  are  observed.   If  the  additional 
normal  stress  be  wanting  in  the  vertical  direction,  the 
boundary  will  tend  to  move  up  or  down  according  to  the 
sign  of  K^  .  The  physical  effect  of  the  additional  normal 
stress  in  the  direction  of  fluid  motion  is  not  as  easily 
interpreted.   We  shall  discuss  this  further  below. 


Equations  of  motion.  -  Recalling  that  fg  (y~ J  *l(*jt)  ,    the 
equations  of  motion  take  the  form 
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6  a/-      **  •       §7      a  a*  vy;  ** ' 


Solutions  by  Cases . 

Case  I .  -  Classical  case.   Suppose  ^J  s%*0,  //  *  £/*' 
and  f    is  constant .   Then  (9-1)  reduces  to 


(9-2) 


O 


M   iS   —  ■  "   • 


Now   from    (9-2) 2»   Jb    is   a  function  of    X    only.      From    (9-1)^ 

to  at  OL  t  X  "S*  CLj 

2 
where  CK»»  and  <Xj,  are  constants.   Boussinesq  studied  the  case 

where  Olj-  O  *  and  pave,  as  a  solution  to  the  flow  equation, 


where 


"'•"  -fl?* 


Case  IIA.  -  The  relationship  between  assumptions  con- 
cerning the  pressure  and  existence  of  flows  of  non-Newtonian 
fluids:,  as  discussed  previously,  is  demonstrated  in  the  prob- 
lem of  Boussinesq  when  3j,  is  considered  constant,  in  which 
case  (9-1)  becomes 

1Ibid. 
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ax* 

<fc-4) 


where,  from  (9-4) 2 

Supposing,  now,  that  we  assume  that  -P   i»  *  constant.  Then 
fro*  the  explicit  relationship  for  £> ,  since  (-ij-r)  is  a 
function  of  Z  and  t  only,  It  must  be  the  case  that 


hence 


(9-6) 


that  is,  the  velocity  is  a  linear  function  of  the  distance 
froa  the  boundary  and  therefore  the  second  derivative  van- 
ishes. The  right  hand  side  of  (9-4 >^  therefore  vanishes,  and 

whence 

U.-kt+/UM  <9-7> 


4  comparison  of  (9-6)  with  the  above  yields 

ft%(t)  *c*ntt.*K 
A*<z)  «  Kz. 


•    ■■  >. 
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so  that  the  velocity  is  of  the  for* 

u,(z,t)~  kt+fa-  <9-7) 

The  boundary  conditions,  applied  to  this  form,  restrict  the 
velocity  field  to  the  trivial  case,  si  identically  aero. 

We  have  therefore  shown  that,  under  assumption  of 
constant  pressure  and  /*  constant,  there  can  be  no  Boussinesq 
flow  of  a  non-Newtonian  fluid. 

Viguier3  discussed  the  problem  of  Boussinesq  in  the 
case  where  the  velocity  gradients  are  not  negligible,  em- 
ploying, in  his  equations  of  motion,  viscosity  forces  as  odd 
functions  of  the  velocity  gradients  M  earlier  introduced 

by  himself. 

Equation  (9-7)  comes  about  as  a  direct  consequence  of 
assuming  that  fr   and  f   are  constant.  Since  *&>**  «0  ,  it  must 
be  the  case,  as  noted  from  the  conclusion  above  regarding  the 
impossibility  of  the  flow,  that  U.{tti)* it   .  Kow  the  assumption 
of  adherence  to  the  boundary  gives  *,fr,t)*6   from  which  must  be 
concluded  that  the  velocity  is  identically  aero.   X  have  inter- 
preted this  result  as  indicating  that  the  assumption  of  the 
constancy  of  pressure  should  not  be  made  and  remove  this 
condition  in  the  cases  studied  below. 


3G.  Viguier,  "Sose  thoughts  on  a  problem  of  M.J. 
Boussinesq,"  Acadjmie  roy^le  de  Belaique  PulleUpn  Cjasse 
des  sciences,  1950-1,  71-76. 

4G.  Viguier,  "The  equation  of  the  boundary  layer  in 
the  case  of  higher  order  velocity  gradients,"  academic  des 
Sciences  Comptes  Rendus,  CCXXIV,  (1947),  713. 
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Vlguier  retains  the  assumption  of  constancy  of  pMMMMrv 
and  removes  the  condition  of  adherence,  stating  "the  condition 
of  adherence  to  the  boundary  does  not  seem  to  be  satisfied* 
However  one  may  attempt  to  give  partial  justification." 

Case  IIB.  -  If  we  assume  that  the  pressure  is  a  function 
of  the  flow  direction  and  time  only,  similar  considerations 
yield  precisely  the  same  results  as  in  Case  IIA  above. 

Case  I1IA.  -  We  are  inevitably  led  in  any  study  of  thf 
problem  of  Boussinesq,  from  the  point  of  view  of  non-Newton 1* it 
fluids,  to  assumption  of  the  dependence  of  the  pressure  ciiwtrx- 
bution  on  the  distance  from  the  boundary.   We  consider  fy     ausd 
£  constant  and  assume  jo  »  jk>(*,t)» 

Then  from  (9-1),  we  obtain, 


t»-i(^r**w 


We  may  take,  as  a  solution  for  the  first  equation,  that  giver* 
by  Bousslnesq 


where,  as  before 

Differentiating  (9-9)  with  respect  to  X  yields 

which  can  be  rewritten,  after  tsttfl-gratiort  by  parts, 


(®..m 
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from  which  the  pressure  may  be  displayed, 

We  consider  the  growth  properties  of  the  pressure.   It 
is  clear  that  the  pressure  is  constant  downstream.   at  the 
boundary , 


f<°'t)=  *Tj?~  +  A*™- 


Returning  to  (9-10) , 
--till***)        „ 

and  the  other  term  in  the  square  bracket  is  easily  seen  to 
vanish  as  £— >  O   .   Hence  M**°)   approaches  n^\°)     as  one  would 
expect.   An  additional  initial  condition  is  necessary  to  establish 
the  function  A±(£)   uniquely. 

Now,  as  Z.  increases  without  bound 

and  therefore  p(Zjt)  approaches  f  (£/   . 

Finally,  as  Z>     increases,  since  J»^  a4jJ  m    jl 

Case  IIIB.  -  Supposing  now  that  the  pressure  is  a  func- 
tion of  X  »  Z  »  and  "t   .   Then 


(9-11) 
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and    ^-^i(^f  ♦#!«,*». 

From  (9-11)  ■£,  since  U    is  independent  of  X  >  -jrr  cannot  be 
a  function  of  X  .   Therefore,  from  (9-11)2* 

and  therefore 

<{>(x,i)  a  xA,(t)f  A*(t). 

Hence  we  obtain,  from  (9-11), 

with  the  integral 

Now  since  (LiO/t^O     and  tne  term  in  the  bracket  vanishes 
at  Z»0»  it  must  be  the  case  that  \  /}i<&t   vanishes  identi- 
cally and  therefore  Ai  vanishes  identically. 

We  point  out ,  finally  that  in  the  completely  general 
case  when  ^^  *  gOi-^elfj  >    *rom  consideration  of  the  general 
explicit  form  of  the  pressure,  we  see  once  more  that  the 
pressure  must  be  a  function  of  the  distance  from  the  boundary 

In  conclusion,  we  have  demonstrated  that  the  second 
order  theory  of  the  flow  of  an  incompressible  non -Newtonian 
fluid  demands,  in  the  problem  of  Boussinesq,  that  the 


t 
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pressure  distribution  must  be  of  specific  functional  de- 
pendence.  It  cannot  be  considered  constant  without  forcing 
the  classical  case  or  a  trivial  case  and,  further,  the  pres- 
sure cannot  be  a  function  of  the  flow  direction. 


X.  THE  DRAINING  OF  A  VERTICAL  PLATE 

General.  -  We  solve  the  equations  of  flow  of  a  flat  plate 
of  infinite  width  which  is  drawn  from  an  infinite  fluid  and 
compute  the  thickness  of  the  fluid  film  on  the  surfaces  of 
the  plate.   We  assume  a  condition  of  steady  state  has  been 
reached,  that  the  fluid  is  incompressible  and  that  the  vel- 
ocity is  vertical  and  independent  of  the  height. 

This  problem  was  first  considered  by  Goucher  and 

1  2 

Ward  and  an  equivalent  problem  studied  by  Jeffreys  . 

We  consider  a  rectangular  cartesian  coordinate  system 
with  the  Uf-axls  horizontal  and  perpendicular  to  the  vert- 
ical plate]  the  Z.  -axis  is  vertical  and  parallel  to  the 
plate,  so  that  when  M»T  ,  the  thickness  of  the  film,  uf*0 
and  ux'*0;   and  that  aT*  V  .  a  constant  (i.e.,  the  vertical 
velocity  of  the  plate),  when  H*0  ,  uT  being  a  function  of 
Of   only.   We  disregard  the  curvature  of  the  surface  in  the 
4t~£  plane  and  assume  the  acceleration  of  the  fluid  is 
negligible. 

Under  the  above  assumptions,  the  velocity  vector  is 
given  by, 

F.S.  Goucher  and  H.  Ward,"  A  problem  in  viscosity: 
the  thickness  of  liquid  films  on  solid  surfaces  under  dynamic 
conditions,"  Philosophical  Magazine  and  Journal  of  Science, 
XLIV,  (1922),  239-264  &  X002-1014. 

2 

H.  Jeffreys,  "The  draining  of  a  vertical  plate," 

Proceedings  of  the  Cambridge  Philosophical  Society,  XXVI, 

(1930)  3  204-205 z=x 

This  is  equivalent  to  assuming  that  there  are  no 
normal  stresses  on  the  free  surface. 
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The  boundary  coicdttiBons  are:: 

ur'(r)  ■■  o 
or  ^o)  **  if 

the  external   fosio*  Jfiild   is 

and  the  rate     of  deformation  tensor   is 

&      o       o 


c5Ui|)    *» 


0        O'     i  <*' 


The  stress  tena*w   takes  tHe  form 

t(ij)  as  -Jo  *  ii  o 

]p    O    I J 


To  o   ©1 
h  i*ar'  o  ©  i  I 

n   ,,•  Toi  ©i  ol     .     ,».,   To  O  ©1 


where 


and  the  Poynting  effect  is  noted.   An  additional  stress  along 
the  ck -axis  is  natedj;  if  this  stress  be  wanting,  the  thickness 
of  the  film  wil  J|  tend,  to  increase. 


iM  Equations  of  mofiionu  -  The  equations  of  motion  take  the  form 
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From  (10-1) 2 »  if  4)  Is  constant,  then  ^  flf*  is  constant. 
But  fA  UJ*      is  a  power  series  In  Af  ,  and  therefore  It  mus  t 
be  the  case  that  UT      Is  constant.   Since  Hf    vanishes  when 
u.*T  i  UT/  must  be  identically  zero,  hence  ST  must  be  const *U 
Since  UF-O   when  a* 7*  ,  <e"  must  also  vanish  identically  and  the 
flow  is  trivial.   We  have  demonstrated  that,  in  the  drain  Las; 
of  a  plate  being  removed  from  a  non-Newtonian  fluid,  the  pres- 
sure cannot  be  assumed  constant,  as  in  the  classical  case. 

The  solution.  -  Since  the  pressure  is  Independent  of  X\  w« 
have  from  (10-1) 2 

Integrating  (10-1)3  with  respect  to  X  »  we  obtain 

Comparing  (10-2)  and  (10-3) ,  we  may  write 

f  *  Mr  £§)*.  +:£-/*&/'*  *****£ 

and  ;i 

the  choice  of  At  following  from  the  absence  of  normal  stresses. 
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By  integration  and,  from  the  boundary  condition  for 
US1    ,  we  obtain 

-Liar's  4i(^-T).  <io-4) 

As  in  chapter  V,  we  revert  (10-4) ,  obtaining 

"' -  i&W  **$&  <-*-jU 

which,  upon  integrating  and  substituting  (3-5)2 # 
When   Ml«-0 

«rW«^*  3Z^"T  * <10-6> 

We  may  therefore  combine  (10-5)  and  (10-6)  and  write 

The  first  approximation, 

is  the  solution  of  Goucher  and  Ward. 

We  now  seek  numerical  values  for  the  thickness ,  "7"  » 
in  the  second  approximation,  in  order  to  compare  these  values 
with  those  reported  experimentally  by  Goucher  and  Ward.  To 
this  end,  we  take  only  the  first  two  terms  on  the  right  in 
the  expansion  given  by  (10-6) ,  and  may  therefore  write 
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where    36  3Gm,A8>     and   is  a  constant  dependent  only  on  the 
fluid  much  the  same  as  M  .      We  see    that  <Lt*0     implies  that 
$  *0    and 

'     ~    A, 

The  thickness  in  the  second  approximation.  -  From  (10-7) , 
we  write  the  expression  for  the  square  of  the  thickness 

In  order  to  compare  the  experimentally  measured  thicknesses 
with  those  given  in  (10-8) ,  we  take 


JJL^.ZS   ,  At   -  1+10 


in  c.g.s.  units,  which  fall  within  the  range  of  the  measure- 
ments of  Goucher  and  Ward. 

Using  a  point  from  the  experimental  results,  we  then 
compute  an  value  tor  JB    in  c.g.s.  units  (dimension  -j —  ). 
<j£  is  negative  and  in  the  interval 

-I0  4  §?  *  10*  £  -J 

with  a  value  of  —  8.  d>  K  '0    yielding  a  curve  which  compares 


If  ^  "/^^  instead  of  {("£***)     »  $   wll3L  be  P081*1 
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very  well  with  observation.   We  rewrite  -  <§>  for  $  in  (10-8) 


From  knowledge  of  the  dependence  of  the  thickness ,  I   ,  on  the 
velocity,  y  ,  in  the  classical  case,  »ud  Iron  the  available 
experimental  results ,  we  know  that  the  thickness  increases 
with  an  increase  in  the  velocity.   Therefore  the  expression 
in  (10-9)  with  the  negative  sign  is  that  which  yields  the 
solution  to  the  physical  problem. 

The  following  table  shows  comparative  values  in  the 
classical,  non-linear  and  experimental  cases: 


VELOCITY 

THICKNESS 

linear 

non-linear 

experimental 

0,00 

,0000 

.0000 

.0000 

0.16 

,0075 

.0076 

.0078 

0.64 

0150 

.0151 

.0165 

1.44 

.0225 

,0229 

.0255 

2.56 

0300 

,0311 

.0345 

4.00 

0375 

.0398 

.0438 

5.76 

0450 

.0494 

.0546 

7.84 

0525 

.0610 

.0669 

10.24 

0600 

,0820 

.0822 

The  above  values  of  T   are  plotted  as  a  function  of  fyT*  in 
figure  2. 

How,  in  the  classical  case,  the  thickness  increases 
without  bound,  as  the  velocity  increases.   Intuitively,  we 
would  expect  that,  for  a  given  fluid,  I.e. f  for  a  given 
<jj ,  P    ,  and  14    ,   there  probably  is  a  velocity  beyond  which 
the  assumptions  regarding  laminar  flow  and  boundary  condi- 
tions at  M*X  are  no  longer  valid.   This  is  evident  in  the 
non-linear  theory. 
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TMe  expression  in  (10-9)  no  longer  has  physical  meaning 

itkeni 

V  >  -§£"  do-io) 

tft   Dike  point   where    the   velocity   is   equal   to  the  expression 
m  tt&e  right    in    (10-10)  ,   we  obtain 

iflkeir.1©   \fw,i   is   the  maximura   value  that    can  be  taken   on   by     Y  ,    so 
tthiit:  the  maxiraum   value   attainable   for   the   thickness    is    given 
jly   for  &   particular   fluid 
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APPENDIX 

DERIVATION  OF  THE  EQUATIONS  OF  HYDRODYNAMICS1 

1.   Equation  of  continuity.  -  The  equation  of  continuity  ex- 
presses the  conservation  of  mass.   Suppose  that  a  sufficiently 
smooth  surface,  S  >  includes  a  volume  of  fluid  mass,  V   . 
We  define  the  mass,  rtt  , 

where  J  indicates  integration  over  the  volume  V  ;  £    is  the 
fluid  density  and  ^V  an  element  of  volume.   We  wish  to  ex- 
press the  principle  that  mass  is  neither  created  nor  destroyed* 
Hence  the  time  rate  of  flow  of  mass  from  a  given  stationary 
volume  must  vanish,  except  for  any  mass  which  is  added  or  sub- 
tracted from  the  volume. 

We  let  -yr-  be  the  rate  of  increase  of  mass.  The 
local  mass  flux  across  an  element  of  area  (directed  outwardly) 
oi$%      P©r  unit  time  is  given  b\r  flk*$&i     where  X*  is  the  velo- 
city vector.   Therefore  to  express  conservation  of  mass,  we 

S. 

write, 


4^- *~  k  e*;&*i 

*t    Ms 

where  $   indicates  integration  over  the  surface  5  .   Heace 

Assuming  ^V*  and  g    continuously  different iable,  er 
applying  Green's  Theorem,  we  obtain 


Wff-*-  te**hi)JLV-o 


Following  the  manner  of  Trueedell  in  lectures  on 
Hydrodynamics  at  Indiana  University,  1st  semester,  1951-1952 
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whence  we  obtain  the  Eulerian  Equation  of  Continuity, 

2.   Cauchy's  Laws  of  Motion   (1821,  1828)   Prom  an  axiomatic 
approach  to  continuum  mechanics,  we  have  that  the  force  and 
moment  acting  on  a  body  consists  of  three  parts: 

i.  A  force  and  moment  equal  to  resultant  force  and 
moment,  respectively,  of  a  certain  vector  field,  *t (n)  »  de- 
fined on  the  boundary:,  <5  ,  of  a  body;  thus 


As 


where  £»'**s^  "  *ji>^  are  an  even  permutation  of  /,£,$ 
or:  — /  if  an  odd  permutation  or  O  if  if}>k  are  not  all  dif- 
ferent. We  call  u^y\  the  stress  vector  and  shall  see  that 
dlfcifeiension    2*,*        is 


M 


ii.   *  force  and  moment  equal  to  the  resultant  force 
sjnd  momentum  respect  lively  of  a  certain  density, 

iihlere  P  ,  Q.  are  two  points.   Then 

irtotere  t  S  rf2  ilwl£<2)  over  *ne  whole  space  represents 
1jhe  effect  of  ^  on  a  point  P  as  Q,  passes  through  all 
points  of  the  space. 

iii.   A  force  and  moment  equal  to  the  resultant  force 


. 


1 


. 
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and  Moment  respectively  of  a  certain  moment  density  j*  (**,  X\*t) 
plus  a  certain  moment  density  Xc.(**,Xv*  t )    ;  then 

^  to, 

and  J**?*  +4m  is  *k°  extraneous  force  field. 

We  therefore  have  the  axioms  of  continuum  mechanics 

where  /p*  and  ^^*  are  the  linear  and  angular  momentum 
respectively.   The  subscript  In)    refers  to  the  side  of  a 
particular  surface,  $  ,  on  which  ^^acts. 

3-   Cauchy's  Laws  of  Motion. 

i*  Cauchy's  First  Theorem  Let  5  be  anv  surfatce 
(assumed  regular  enough  to  allow  application  of  Stokes r  and 
Green*  s  Theorems) „  ^JL#i  and  t(ir)  De  tne  vector  dens  it  lies 
on  the  «$■  side  and  —  side  of  the  surface,  5  »  each  contiin- 
uous  on  its  own  side;  then  if  4  -  X*   is  of  bounded  variation, 


0*>  ■  Wa-) 


Proof:   This  theorem  is  proved  from  the  equation  i 
for  JO   applied  to  a  surface  on  each  side  of  $ 


Taking  as  a  bound  on  the  integral  of  -J--Xv   ,  smnd 
allowing  the  surfaces  to  shrink  to  5  »  we  obta::Ln 

whence  from  the  assumed  continuity, 


■ 


- 


• 
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ii.   Cauchy's  Second  Theorem  Using  the  expression 
for  "%?     we  show  the  existence  of  a  stress  tensor,  defined  by 

where  Hi     is  the  unit  normal  to  a  surface  (directed  outward) . 
The  use  of  the  stress  tensor  allows  us  to  describe  the  stresses 
on  a  surface  without  having  to  speak  of  the  stress  direction. 

iii.   Cauchy's  Laws  of  Motion.  -  ffe  consider  only  the 
continuously  differentiable  case.   From  the  expression  for  4p    > 

How, 

by  Green's  Theorem.   Upon  substituting  &m*  PcLV ,  we  obtain 


it 


V 


Cauchy's  First  Law  of  Motion. 

From  the  expression  for  ft    ,   we  can  similarily  show 
that  £,* .£  i*    y.  0 %l   -  O    »  Cauchy's  Second  Law  of  Motion. 

Now  if  Xi  ,    the  extraneous  moment  of  momentum 
vanishes,  then  "£  *-*T         ,  which  we  henceforth  assume. 

4.   The  Energy  Equatuion.  -  We  define  the  kinetic  energy,  /£  , 

Km±fe*i*iM 


Y6 
We  further  consider  the  internal  energy  if  of  a  body,  and 
assume  that  it  is  a  continuous  function  of  mass.   We  may 
then  consider  an  integral  form  £  =  JeMm    where  €.    is  the  in- 
ternal energy  density.  Iz    is  changed  by  two  specific  pro- 
cesses, (1)  mechanical  work  and  (2)  non-mechanical  causes 
(e.g.  radiation,  heat  conduction,  etc.) 

Principle  of  Energy.  -  Let  there  be  a  heat  flux, 
#*  ,  (analogue  of  the  stress  vector)  such  that  we  may  write 

that  is,  the  total  time  rate  of  change  of  energy  is  equal 
to  the  mechanical  work  performed  plus  the  work  done  due  to 
the  stress  tensor.  <b  «*<$$£   is  so  defined  as  to  balance  the 
equation.   We  note  that  the  interconvertability  of  heat  and 
work  is  already  determined  by  this  equation  because  of  the 
fact  that  Jf*  is  written  in  mechanical  units,  i.e. 

dimension  :a.   energy    m     M 
time  area      TTJ 

Assuming  all  functions  involved  are  continuously  differenti- 
able,  substituting  the  expression  for  |^  and  £  and  applying 
the  divergence  theorem,  we  obtain 

Applying  CJauchy's _ first  law,  and  from  continuity  assumptions, 


we  have 


•  * 


f£«r'  *i.i-%9*' 
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Defining  the  rate  of  deformation  tensor  as  usual, 

J 


aij^i-ix^-t-Kj.i) 


and  recalling  the  symmetry  of  the  stress  tensor,  we  have 

«         3    *  t 

P£*tCicl*«*   -  9*>i        "the  Fourier -Kirchof f -C.   Neumann  equation 
Thus  £lM<lg   ,    the  stress  power,  is  the  rate  at  which  stress 
does  work  per  unit  volume. 
Assuming  as  usual  for  a  homogeneous  fluid  that  E    is  given 


by  <£  *  £(V^mJ,  where  Vs  y  and  #  is  a.  local  statistical 

parameter  (entropy),  we  define  7?s  —jjjj  1 0 s  jr  >  the  pressure 

and  temperature  respectively;  then  iA£  =  -7r^V  ■+&&*}. 

We  now  define  the  extra  stress  tensor  V"j  by 
« 

Vj  »  £^1  *  t  i  »  where  >p  is  any  scalar.   After  some  al- 
gebraic manipulation,  we  obtain 


Hence 


gi  -  -f>Jtk  *  r**d*j  -  f*,£ 


&  fj  =  £  rrir  -h  £$  +  Y**Jx  ~  9  \l 


From  the  equation  of  continuity 


whence 

For  an  incompressible  fluid  y"  is  O  ,  and 

*  *    ft  * 
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